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(57) Abstract 

An interpolation circuit capable of performing interpolation operation with a simple constitution. A I6-times overeampling from 
discrete data is performed by D flip-flops (4, 5). A first convolution operation is performed by D flip-flops (4-11) and an adder (12) on the 
result of which a second convolution operation is performed by D flip-flops (13-20) and an adder (21). Data interpolated along a quadratic 
function curve enveloping the discrete data is obtained from the adder (21). 
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DESCRIPTION 



INTERPOLATION CIRCUIT 

TECHNICAL FIELD 

The present Invention relates to an interpolation circuit 
for generating interpolated data from discrete data. 

BACKGROUND ART 

Conventional CD (Compact Disk) reproducing devices 
employ a digital-analog converter of the oversampling method. 
This digital-analog converter performs interpolation 
operation for digital data discretely input, employing a sine 
function or the like. However, this sine function converges 

to 0 at ±oo f and has a shortcoming that if the arithmetic 
operation is truncated at the finite value, a truncation error 
occurs in the operation. Generally, a stepwise waveform 
obtained in the interpolation operation is passed through a 
low-pass filter, but there is also a drawback that a phase 
distortion or a distorted output waveform may arise through 
the low-pass filter. 

Particularly, in the above interpolation operation using 
the sine function, a method is taken in which the values of 
the sine function are stored in a table, and read as needed, 
or held as the tap coefficients for a digital filter, but this 
method has a problem that the configuration is complex. 




Therefore, an interpolation circuit capable of performing 
interpolation operation with a simple constitution is desired. 

DISCLOSURE OF THE INVENTION 

The present invention has been achieved in the light of 
the aforementioned problems, and it is an object of the 
invention to provide an interpolation circuit capable of 
performing interpolation operation with a simple 
constitution. 

The interpolation circuit of the invention comprises 
oversampling operation unit for performing oversampling 
operation from the zero-order hold input data, and first 
convolution operation unit for performing convolution 
operation on plural data obtained by the oversampling operation 
unit twice or more repeatedly, characterized in that 
interpolated data is generated along a quadratic function curve 
passing through the integral multiples of the value of input 
data. 

Also, the interpolation circuit of the invention 
comprises oversampling operation unit for performing 
oversampling operation from the zero -order hold input data, 
second convolution operation unit for performing convolution 
operation on plural first data obtained by the oversampling 
operation unit to calculate plural second data enveloped by 
a symmetrical trapezoid of the shape having an upside of 
substantially 0 . 5 times the width of original input data and 
a base of substantially 1.5 times the width, and third 
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convolution operation unit for performing convolution 
operation on the plural second data obtained by the second 
convolution operation unit to calculate plural third data 
enveloped by a smooth quadratic function curve having a width 
of the base being substantially twice that of the original 
input data. 

Also , the interpolation circuit of the invention 
comprises oversampling operation unit for performing 
oversampling operation from the zero-order hold input data, 
and fourth convolution operation unit for performing 
convolution operation on plural data obtained by the 
oversampling operation unit to calculate plural data enveloped 
by an equilateral triangle having a base being substantially 
twice the width of original input data. 

Also, the interpolation circuit of the invention 
comprises oversampling operation unit for performing 
oversampling operation from the zero -order hold input data 
with a sampling period of 2nTi at a time interval of Ti, fifth 
convolution operation unit for performing convolution 
operation of n phases by adding plural data obtained by the 
oversampling operation unit n times with the data shifted by 
the time interval of T lt and sixth convolution operation unit 
for performing convolution operation of n phases by adding 
plural data obtained by the fifth convolution operation unit 
n times with the data shifted by the time interval of Ti. 

In particular, at least one of the fifth and sixth 
convolution operation unit desirably comprises data holding 



unit for holding n pieces of data output from the oversampling 
operation unit while shifting, and addition unit for adding 
n pieces of data held in the data holding unit. 

In the former stage of the oversampling operation unit, 
it is desirable to provide data appending unit for appending 

the data having the symmetrical values proportional to the 

<} 

input data before and after the input data. 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 is a graph for explaining a principle of 
interpolation operation that is performed by digital -analog 
converter of the first embodiment; 

FIG. 2 is a graph for explaining a principle of 
interpolation operation that is performed by digital -analog 
converter of the first embodiment; 

FIG. 3 is a graph for explaining a principle of 
interpolation operation that is performed by digital -analog 
converter of the first embodiment; 

FIG. 4 is a chart showing a specific example of an 
interpolation operation; 

FIG. 5 is a chart showing a specific example of an 
interpolation operation ; 

FIG. 6 is a graph showing a waveform that is finally 
obtained by arithmetic operations as shown in FIGS. 4 and 5; 

FIG. 7 is a diagram showing a waveform corresponding to 
the partial data shown in FIG. 4; 
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FIG. 8 is a diagram showing a process for performing the 
first convolution operation; 

FIG. 9 is a diagram showing a process for performing the 
second convolution operation; 

FIG. 10 is a diagram showing a process for performing 
the third convolution operation; 

FIG. 11 is a diagram showing a specific example of the 
interpolation operation when 0, 3, 7, 5, -4, and 0 are input 
as the discrete data; 

FIG. 12 is a graph showing a waveform that is obtained 
by arithmetic operation as shown in FIG. 11; 

FIG. 13 is a diagram showing the configuration of a 
digital- analog converter of the first embodiment; 

FIG. 14 is a diagram showing the configuration of a 
digital-analog converter of the second embodiment; 

FIG. 15 is a diagram showing a specific example of 
arithmetic operation in the first , second and third convolution 
operation circuits included in the digital- analog converter 
as shown in FIG. 14; 

FIG. 16 is a graph showing a waveform that is finally 
obtained by arithmetic operation as shown in FIG. 15; 

FIG. 17 is a diagram showing a waveform with another pulse 
symmetrically appended on both sides of the single input pulse; 

FIG. 18 is a graph showing a result of effecting plural 
convolution, after appending another pulse of FIG. 17 to the 
input data; 
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FIG. 19 is a circuit diagram showing the configuration 
of a circuit for appending another pulse on both sides of the 
zero-order hold, single pulse; 

FIG. 20 is a graph for explaining the principle of linear 
interpolation ; 

FIG. 21 is a circuit diagram showing the configuration 
of a digital -analog converter for linear interpolation of the 
fifth embodiment; 

FIG. 22 is a diagram showing a specific example of 
interpolation operation in the case where the discrete input 
data are 0, 3, 7, 5, -4 and 0; and 

FIG. 23 is a graph showing a waveform that is obtained 
by arithmetic operation as shown in FIG. 22. 

BEST MODE FOR CARRYING OUT THE INVENTION 

Embodiments of a digital-analog converter to which a data 
interpolation method of the present invention is applied will 
be described below with reference to the accompanying drawings . 

[First embodiment] 

A digital-analog converter in a first embodiment of the 
invention performs interpolation operation employing a 
quadratic function curve connecting the input data smoothly 
when the discrete digital data is input. 

FIGS. 1 through 3 are graphs for explaining the principle 
of interpolation operation that is performed by the 
digital -analog converter in this embodiment of the invention, 



in which the single data is converted into data interpolated 
along a quadratic function curve. As shown in FIG. 1, the 
zero-order hold single data is considered to have an amplitude 
of 1 and a length of 2n-Ti. When the discrete data is input, 
the length 2n-Ti for the single data is a cycle period (sampling 
frequency) of input data. For the data as shown in FIG. 1, 
the first convolution operation is performed by shifting the 
data by T L and performing oversampling operation from the data , 
and adding the data n times. Consequently, a symmetrical 
trapezoidal wave with the base of (3n-l)*T 1# the upside of 
(n+l)*Ti and the height of n can be obtained, as shown in FIG. 
2. Further, for this symmetrical trapezoidal wave, a second 
convolution operation is performed by shifting the data by 
Ti and adding the data n times, so that a continuous quadratic 
function curve with the width of (4n-l)*Ti and the amplitude 
of 2n 2 can be obtained, as shown in FIG. 3. 

In this way, by performing oversampling operation for 
the zero- order hold input data and then performing convolution 
operation of n phases , the output data enveloped by the 
continuous quadratic function curve can be obtained. 
Therefore, considering a case of zero -order holding the 
discrete input data successively, the quadratic function curve 
corresponding to each input data is generated, while being 
shifted by 2nTi. As a result, the output is 2n pieces of data 
interpolated along the quadratic function curve passing 
through and smoothly connecting 2n 2 times the input data. 



FIGS. 4 and 5 are charts showing a specific example of 
the above mentioned interpolation operation. Herein, a 
calculation example with n=8 is shown. The detailed contents 
of the first convolution operation are shown in FIG. 4, and 
the detailed contents of the second convolution operation are 
shown in FIG. 5. 

In FIG. 4, each item of data listed at (1) is the zero-order 
hold data (data after oversampling operation) as shown in FIG. 
1 . The elapse of time is arranged in the transversal direction , 
the unit of time being Ti corresponding to one column. Namely, 
it is supposed that the input data corresponding to FIG. 1 
are 16 pieces of zero-order hold data, with the duration of 
16Ti and an amplitude of 1 . An 8-phase convolution operation 
is performed, comprising producing eight sets of data by 
shifting the zero-order hold data by time T i# as shown in (1) 
through (8) in FIG. 4, and adding them. Consequently, a result 
of calculation is obtained as listed at (Ti) in FIG. 4. From 
this result of calculation, it will be found that a symmetrical 
trapezoidal wave as shown in FIG. 2 is obtained. 

In FIG. 5, each item of data listed at (1) is the data 
obtained by the first convolution operation, and the same as 
shown at (Ti) in FIG. 4. As listed at (1) through (8) in FIG. 
5, an 8-phase convolution operation is performed, comprising 
producing eight sets of data by shifting the data by time Ti, 
as shown in (1) through (8) in FIG. 5, and adding them. 
Consequently, a result of calculation is obtained as listed 
at (T2) in FIG. 5. (Al) of FIG. 5 is the data obtained by 
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the second convolution operation, and the same as listed at 
(T2). A two-phase convolution operation (addition of (Al) 
data and (A2) data) is performed, comprising shifting the data 
listed at (Al) by T x ((A2) in FIG. 5), so that a result of 
calculation is obtained as listed at (A3) in FIG. 5. 

FIG. 6 is a graph showing a waveform that is finally- 
obtained by arithmetic operations as shown in FIGS. 4 and 5. 
The longitudinal axis indicates the amplitude, and the 
transversal axis indicates the time, respectively. As shown 
in FIG. 6, an original rectangular wave of zero-order hold 
data results in the oversampled data enveloped by a smooth 
quadratic function curve with a double width, the maximum 
amplitude being 128 times (= 2n 2 ) . 

The above example is described employing the single input 
data. However, in the case where plural pieces of data are 
successively input , the data oversampled corresponding to each 
input data are arranged without overlapping each other. 
Therefore, when the data is input successively at a period 
of 2n-Ti, an interpolation curve always passes through the 
integral multiple (128 times) of each item of data. 

FIG. 7 is a diagram showing a waveform corresponding to 
the data of (1) in FIG. 4. This waveform corresponds to a 
zero-order hold unit pulse with an amplitude of 1 and a width 
of 2n-Ti. FIG. 8 is a diagram showing a process for performing 
the first convolution operation, and schematically showing 
how the amplitude and the width change. Also, FIG. 9 is a 
diagram showing a process for performing the second convolution 




operation, and schematically showing how the amplitude and 
the width change. FIG. 10 is a diagram showing a process for 
performing the third convolution operation , and schematically 
showing how the amplitude and the width change . As shown in 
FIG. 10, the third convolution operation results in the total 
width (data length) of plural data after interpolation being 
(4n-l)-T L . 

Then, a continuity test for the waveform obtained by the 
second convolution operation is made. In FIG. 5, with the 
position in the transversal direction corresponding to x, the 
value y of each item of data at ( A3 ) is obtained in the following 
way. 

For 0 <; x * n, (y in this interval is denoted as yi) 
y x = (1+2+ ... +x)+(l+2+ ... +(x-l)) 
- x(x+l)/2+(x-l)x/2 

= x 2 ... (1) 

For n ^ x £ 3n, (y in this interval is denoted as y 2 ) 
y 2 = 2(1+2+ ... + n )-(l+2+ ... +(x-n)) 

+n(2x-2n-l) 

= n(n+l)-(x-n) 2 +n(2x-2n-l) 

= -x 2 +4nx-2n 2 ... (2) 

For 3n s x ss 4n, (y in this interval is denoted as y 3 ) 

y 3 = (-x+4n) 2 

= x 2 -8nx+16n ... (3) 

In this way, it will be found that y is in a quadratic expression 
of x. 

The inclination at x = n is obtained as follows . 




dyi/dx = 2x = 2n 

dy 2 /dx = -2x+4n = 2n 
Since both values are equal to 2n, indicating the inclinations 
at x = n are the same, it will be found that the data is smoothly 
connected. 

Also, the inclination at x=3n is obtained as follows. 

dy 2 /dx = -2x+4n = -2n 

dy 3 /dx = 2x-8n = -2n 
Since both values are equal to -2n, indicating that the 
inclinations at x=3n are the same, it will be found that the 
data is smoothly connected. 

FIG. 11 is a diagram showing a specific example of the 
interpolation operation when 0, 3, 7, 5, -4, and 0 are input 
as the discrete data. For example, the specific example is 
shown in a case of n=4. 

As listed at (1) in FIG. 11, first of all, eight zero-order 
hold oversampled data are generated from the input data. A 
four-phase convolution operation is performed, comprising 
generating four sets of data by shifting the eight pieces of 
input data by time T i# as listed at (1) through (4) in FIG. 
11, and then adding them, whereby a result of calculation can 
be obtained as listed at (5) in FIG. 11. 

Similarly, a four-phase convolution operation is 
performed, comprising generating four sets of data by shifting 
the first convolution operation result as obtained in this 
way by time T i# as listed at (5) through (8) in FIG. 11, and 
then adding them, whereby a result of calculation can be 



obtained as listed at (Al) in FIG, 11. Then a two-phase 
convolution operation (addition of (Al) data and (A2) data) 
is performed, comprising further shifting the data listed at 
(Al) by Ti ((A2) in FIG. 11). Consequently, a calculation 
result is obtained as listed at (A3) in FIG. 11. 

FIG. 12 is a graph showing a waveform that is obtained 
by arithmetic operation as shown in FIG. 11 . The longitudinal 
axis indicates the amplitude, and the transversal axis 
indicates the time. As shown in FIG. 12, it can be found that 
the discrete input data is interpolated by the data connect 
with these data by quadratic function curve. The amplitude 
of output value corresponding to each input data is 32 times 
(=2n 2 ) the value of input data. 

FIG. 13 is a diagram showing the configuration of a 
digital -analog converter in this embodiment of the invention. 
For example, a specific configuration is shown in a case of 
n=8 . 

As shown in FIG. 13, the digital-analog converter of this 
embodiment comprises eight D flip-flops 4 through 11 and an 
adder (ADD) 12 making up a first convolution operation circuit , 
eight D flip-flops 13 through 20 and an adder (ADD) 21 making 
up a second convolution operation circuit, a D flip-flop 22 
and an adder (ADD) 23 making up a third convolution operation 
circuit, a divider 24, a digital-analog converter (D/A) 25, 
and a filter 26. 

A basic clock signal CLK1 has the same frequency as the 
sampling frequency of the discrete input data, and is input 




into an initial stage D flip-flop 4, The discrete input data 
is fetched into the D flip-flop 4 in synchronism with this 
basic clock signal CLK1 , and held. Also, a clock signal CLK2 
has a frequency of 16 times the basic clock signal CLK1 . All 
the D flip-flops at the second stage and beyond perform the 
operation of sampling the data in synchronism with the clock 
signal CLK2 . The oversampling operation unit is constituted 
by two D flip-flops 4, 5 into which the clock signals CLK1, 
CLK2 having different f requencies are input . 

If 16 -bit data is input into the initial stage D flip-flop 
4 in the first convolution operation circuit, seven D 
flip-flops 5 through 11 at the second stage and beyond accept 
the data held in this initial stage D flip-flop 4 in synchronism 
with the clock signal CLK2 in order and successively and shift 
the data. The adder 12 adds the data held in the eight D 
flip-flops 4 through 11. In this way, a result of arithmetic 
operation as listed at (Ti) in FIG. 4 can be obtained in the 
adder 12 within the first convolution operation circuit. 

Also, eight D flip-flops 13 through 20 included in the 
second convolution operation circuit accept the 19 -bit data 
output from the adder 12 within the first convolution operation 
circuit in synchronism with the clock signal CLK2 in order 
and successively and shift the data. The adder 21 adds the 
data held in the eight D flip-flops 13 through 20. In this 
way, a result of arithmetic operation as listed at (T2) in 
FIG. 5 can be obtained in the adder 21 within the second 
convolution operation circuit. 



Further, the D flip-flop 22 included in the third 
convolution operation circuit accepts the 22 -bit data output 
from the adder 21 within the second convolution operation 
circuit in synchronism with the clock signal CLK2 . The adder 
23 adds the data output from the adder 21 and the data held 
in the D flip-flop 22. In this way, a result of arithmetic 
operation as listed at (A3) in FIG. 5 can be obtained in the 
adder 23 within the third convolution operation circuit. 

The 2 3 -bit data output from the adder 23 is divided by 
8 in the divider 24, and converted into the 20-bit data. This 
20-bit data is converted into a stepwise analog signal by the 
digital -analog converter 25, which is then passed through the 
filter 26 for suppressing unnecessary radiation by removing 
a frequency component of 16 times the basic clock signal CLK1 , 
and output . 

[Second embodiment] 

By the way, in the first embodiment as described above, 
2n pieces of zero -order hold data are produced by the 
oversampling operation on the basis of the input data and 
shifted by the data interval of T x , and n sets of shifted data 
are added. However, the circuit can be simplified by modifying 
the addition order. 

FIG. 14 is a diagram showing the configuration of a 
digital -analog converter according to a second embodiment of 
the invention . For example , a specific configuration is shown 
in a case of n=*8 . 



As shown in FIG, 14, the digital -analog converter of this 
embodiment comprises nine D flip-flops 31 through 36 , 38, 39, 
and 41 and three adders (ADD) 37, 40, 42 making up a first 
convolution operation circuit, seven D flip-flops 43 through 
46, 48, 49 and 51 and three adders (ADD) 47, 50, 52 making 
up a second convolution operation circuit, a D flip-flop 53 
and an adder (ADD) 54 making up a third convolution operation 
circuit, a divider 55, a digital -analog converter (D/A) 56, 
and a filter 57. 

Instead of adding eight sets of data that are obtained 
by shifting eight pieces of zero-order hold data by the data 
interval T 1# the first convolution operation circuit involves 
a first addition operation of adding two sets of data shifted 
by the interval 4Ti at a time, a second addition operation 
of adding two sets of data shifted by the interval 2T X at a 
time, and a third addition operation of adding two sets of 
data shifted by the interval Ti . 

More specifically, employing f our D flip-flops 33 through 
36 included in the first convolution operation circuit, the 
data shifted by 4Tx (i.e., four periods of the clock signal 
CLK2) are generated from eight pieces of zero-order hold data 
successively input , and the first addition operation of adding 
the data output from the D flip-flop 32 and the data output 
from the D flip-flop 36 is performed by the adder 37. Also, 
employing two D flip-flops 38, 39, the data shifted by 2Ti 
are generated from the data output from the adder 37, and the 
second addition operation of adding the data output from the 



adder 37 and the data output from the D flip-flop 39 is performed 
by the adder 40. Further, employing the D flip-flop 41, the 
data shifted by Ti is generated from the data output from the 
adder 40, and the third addition operation of adding the data 
output from the adder 40 and the data output from the D flip-flop 
41 is performed by the adder 42. In this way, the same 
convolution operation is performed in the same manner as the 
first convolution operation circuit constituted by the eight 
D flip-flops 4 through 11 and the adder 12 as shown in FIG. 
13. 

Similarly, employing four D flip-flops 43 through 46 
included in the second convolution operation circuit, the data 
shifted by 4Ti are generated from the data output from the 
adder 42 within the first convolution operation circuit, and 
the first addition operation of adding the data output from 
the adder 42 and the data output from the D flip-flop 46 is 
performed by the adder 47. Also, employing two D flip-flops 
48, 49, the data shifted by 2Ti is generated from the data 
output from the adder 47, and the second addition operation 
of adding the data output from the adder 47 and the data output 
from the D flip-flop 49 is performed by the adder 50 . Further, 
employing the D flip-flop 51 , the data shifted by Ti is generated 
from the data output from the adder 50, and the third addition 
operation of adding the data output from the adder 50 and the 
data output from the D flip-flop 51 is performed by the adder 
52. In this way, the same convolution operation is performed 
in the same manner as the second convolution operation circuit 
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constituted by the eight D flip-flops 13 through 20 and the 
adder 21 as shown in FIG. 13. 

Further, the D flip-flop 53 included in the third 
convolution operation circuit accepts the data output from 
the adder 52 within the second convolution operation circuit 
in synchronism with the clock signal CLK2. The adder 54 adds 
the data output from the adder 52 and the data output from 
the D flip-flop 53. In this way, the convolution operation 
is performed by the third convolution operation circuit . 

The 23 -bit data output from the adder 54 is divided by 
8 in the divider 55, and converted into the 20 -bit data. This 
20 -bit data is then converted into a stepwise analog signal 
by the digital-analog converter 56, which is then passed 
through a filter 57 for suppressing unnecessary radiation by 
removing a frequency component of 16 times the basic clock 
signal CLK1, and output. 

FIG. 15 is a diagram showing a specific example of 
arithmetic operation in the first , second and third convolution 
operation circuits included in the digital -analog converter 
as shown in FIG. 14. 

In FIG. 15, each item of data as listed at (1) shows the 
input data from the D flip-flop 32 into the adder 37, and each 
item of data as listed at (2) shows the data input from the 
D flip-flop 36 into the adder 37. Also, each item of data 
as listed at (Ti) and (3) is the same, and is output from the 
adder 37. Each data as listed at (4) shows the data input 
from the D flip-flop 39 into the adder 40. Also, each item 



of data as listed at (T2) and (5) is the same, and is output 
from the adder 40. Each data as listed at ( 6 ) shows the data 
input from the D flip-flop 41 into the adder 42. Also, each 
item of data as listed at (T3) and (7) is the same, and is 
output from the adder 42. In this way, it will be found that 
the data ( (T3) in FIG. 15) obtained by the first convolution 
operation circuit is exactly the same as the result of 
convolution operation as listed at (Ti) in FIG. 4. 

Similarly, in FIG. 15, each item of data as listed at 
(7) shows the data input from the adder 42 into the adder 47, 
and each item of data as listed at ( 8 ) shows the data input 
from the D flip-flop 46 into the adder 47. Also, each item 
of data as listed at (T4) and (9) is the same, and is output 
from the adder 47. Each data as listed at (10) shows the data 
input from the D flip-flop 49 into the adder 50. Also, each 
item of data as listed at (T5) and (11) is the same, and is 
output from the adder 50. Each data as listed at (12) shows 
the data input from the D flip-flop 51 into the adder 52 . Also, 
each item of data as listed at (T6) and (Al) is the same, and 
is output from the adder 52. In this way, it will be found 
that the data ((T6) in FIG. 15) obtained by the second 
convolution operation circuit is exactly the same as the result 
of convolution operation as listed at (T2) in FIG. 5. 

FIG. 16 is a graph showing a waveform that is finally 
obtained by arithmetic operation as shown in FIG. 15. The 
longitudinal axis indicates the amplitude , and the transversal 
axis indicates the time, respectively. As shown in FIG. 16, 




an original rectangular wave of zero-order hold data results 
in the oversampled data enveloped by a smooth quadratic 
function curve with a double width , the maximum amplitude being 
128 times ( = 2n 2 ) . 

In this way, the configuration can be greatly simplified 
by changing the order of addition. For example, two adders 
12, 21 having eight-input are employed in the configuration 
as shown in FIG. 13, but they can be implemented by seven 
two-input adders, or a total of 14 adders. On the other hand, 
six two-input adders 37, 40, 42, 47, 50 and 52 are employed, 
instead of two eight-input adders 12, 21 in the configuration 
as shown in FIG. 14, and two adders 12, 21 having eight-input 
can be removed. 

[Fourth embodiment] 

By the way, in the above embodiments, the discrete input 
data is considered to be a single pulse as shown in FIG. 1. 
However, a pulse with different amplitude and polarity may 
be symmetrically appended on both sides (fore and back sides) 
of this single pulse. 

FIG. 17 is a diagram showing a waveform with another pulse 
symmetrically appended on both sides of the single input pulse . 

FIG. 18 is a graph showing a result of effecting plural 
convolution operations on the basis of the operational 
principle as shown in FIGS . 1 through 3 , after appending another 
pulse of FIG. 17 to the input data. As shown in FIG. 18, a 
typical interpolation function having positive and negative 
values can be obtained by performing plural convolution 




operations after appending a pulse with different amplitude 
and polarity on both sides of the single pulse. 

FIG . 19 is a diagram showing the configuration of a circuit 
for appending another pulse on both sides of the zero-order 
hold, single pulse . A pulse appending circuit (corresponding 
to data appending unit) as shown in FIG. 19 comprises six D 
flip-flops 60 through 64, 72, four multipliers 65 through 68 , 
and three adders 69 through 71. An initial stage D flip-flop 
60 fetches the data in synchronism with the basic clock signal 
CLK1 and holds it. Four D flip-flops 61 through 64 connected 
in tandem operate in synchronism with a clock signal CLK4 having 
a frequency of one-half the basic clock signal CLK1, to accept 
the data output from the initial stage D flip-flop 60 in order 
and successively and hold it. Each output data of the D 
flip-flop 61, 64 is input into the multiplier 65, 66 having 
a multiplying factor ™-5", respectively. Each output data 
of the D flip-flop 62, 63 is input into the multiplier with 
a multiplying factor "31". And the results of multiplication 
from each of these four multipliers 65 through 68 are added 
by three adders 69 through 71. In this way, a pulse of FIG. 
17 is output from a final stage adder 71. It should be noted 
that various interpolation functions can be produced by 
changing the number of stages for the D flip-flop and the 
multiplying factor or polarity of the multiplier. 

[Fifth embodiment] 

In the first and second embodiments as previously 
described, the discrete input data is interpolated between 
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the data by means of a quadratic function curve. However, 
the data may be interpolated using the straight line for various 
purposes . 

FIG. 20 is a graph for explaining the principle of linear 
interpolation. A convolution operation of shifting the 
zero -order hold data with a length of 2n*Ti and an amplitude 
of 1 as shown in FIG. 1 by Ti and adding them 2n times, so 
that an equilateral triangular wave with a base of (4n~l)-Ti 
and an amplitude of 2n as shown in FIG. 20 can be obtained. 

In this way, the output data enveloped by the equilateral 
triangular wave can be obtained by over sampling the zero -order 
hold input data 2n times and performing the 2n-phase 
convolution operation. Accordingly, considering the case of 
zero-order holding the discrete input data successively, the 
equilateral triangular wave corresponding to each input data 
is generated, shifted by 2n*Ti, and the output data is 2n pieces 
of data passing through 2n times the input data and connecting 
them along the equilateral triangular wave. 

FIG. 21 is a diagram showing the configuration of a 
digital -analog converter for linear interpolation in this 
embodiment. For example, a specific configuration in a case 
of n=4 is shown . 

As shown in FIG. 21, the digital-analog converter of this 
embodiment comprises eight D flip-flops 81 through 88 and an 
adder (ADD) 89 making up a convolution operation circuit, a 
divider 90, a digital-analog (D/A) converter 91, and a filter 
92. 



If the data is input into an initial stage D flip-flop 
81 included in this convolution operation circuit, seven D 
flip-flops 82 through 88 at the second stage and beyond fetch 
the data held in this initial stage D flip-flop 81 in synchronism 
with a clock signal CLK5 (having a frequency of eight times 
the basic clock signal CLK1) in succession and shift them. 
The adder 88 adds the data held in eight D flip-flops 81 through 
88, respectively. In this way, the convolution operation is 
performed. 

The data output from the adder 88 is divided by 8 in the 
divider 89, and then converted into a stepwise analog signal 
by the digital-analog converter 91, which is then passed 
through a filter 9 2 for suppressing unnecessary radiation by 
removing a frequency component of eight times the basic clock 
signal CLK1 , and output . 

FIG. 22 is a diagram showing a specific example of 
interpolation operation in the case where the discrete input 
data are 0, 3, 7, 5, -4 and 0. For example, a specific example 
in a case of n=4 is shown. 

As listed at (1) in FIG. 22, first of all, eight pieces 
of zero-order hold data is generated for the input data. An 
eight -phase convolution operation is performed, comprising 
generating eight sets of input data by shifting the eight pieces 
of input data by time T 1# as listed at (1) through (8) in FIG. 
22, and then adding them in the adder 89, so that a result 
of calculation as listed at (T) in FIG. 22 can be obtained. 
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FIG. 2 3 is a graph showing a waveform that is obtained 
by arithmetic operation as shown in FIG . 22 . The longitudinal 
axis indicates the amplitude, and the transversal axis 
indicates the time. As shown in FIG. 23, it can be found that 
the discrete input data is interpolated by data connecting 
with these data by a straight line. The amplitude of output 
value corresponding to each input data is eight times (=2n) 
the value of input data . The digital-analog converter as shown 
in FIG. 22 restores the amplitude to an original value, 
employing a divider 90. 

In the above embodiments, this invention is applied to 
the digital -analog converter. However, an oversampling 
circuit may be configured in which an oversampling operation 
of 2n times the input data is performed by removing the divider, 
the digital -analog converter and the filter included in each 
digital-analog converter. 

INDUSTRIAL APPLICABILITY 

As described above, with this invention, the zero- order 
hold discrete data can be interpolated with a simple 
constitution to obtain a smooth waveform or a linearly 
interpolated waveform. For example, the digital -analog 
converter of this invention may be replaced with the 
conventional digital-analog converter used for the digital 
audio equipment by adding an oscillation circuit with PLL to 
improve the quality of sound. 




Since the interpolated data can be easily generated 
without referring to the table in the image processing, an 
extremely wide application range can be considered. 
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CAIMS 

1. An interpolation circuit comprising: 
oversampling operation unit for performing oversampling 

operation from zero-order hold input data; and 

first convolution operation unit for performing 

convolution operation on plural data obtained by said 

oversampling operation unit twice or more repeatedly, 

characterized in that interpolated data is generated 

along a quadratic function curve passing through the integral 

multiples of the value of said input data. 

2. An interpolation circuit, characterized by comprising: 
oversampling operation unit for performing oversampling 

operation from zero -order hold input data; 

second convolution operation unit for performing 
convolution operation on plural first data obtained by said 
oversampling operation unit to calculate plural second data 
enveloped by a symmetrical trapezoid of the shape having an 
upside of substantially 0.5 times a width of original input 
data and a base of substantially 1.5 times the width; and 

third convolution operation unit for performing 
convolution operation on said plural second data obtained by 
said second convolution operation unit to calculate plural 
third data enveloped by a smooth quadratic function curve 
having the width of a base being substantially twice the width 
of original input data. 

3. An interpolation circuit, characterized by comprising: 



oversampling operation unit for performing overs ampling 
operation from zero- order hold input data; and 

fourth convolution operation unit for performing 
convolution operation on plural data obtained by said 
oversampling operation unit to calculate plural data enveloped 
by an equilateral triangle having the width of a base being 
substantially twice the width of said input data. 

4. An interpolation circuit, characterized by comprising: 
oversampling operation unit for performing oversampling 

operation from zero-order hold input data with a sampling 
period of 2n-Ti at a time interval of Ti; 

fifth convolution operation unit for performing 
convolution operation of n phases by adding plural data 
obtained by said oversampling operation unit n times with the 
data shifted by the time interval of T x ; and 

sixth convolution operation unit for performing 
convolution operation of n phases by adding plural data 
obtained by said fifth convolution operation unit n times with 
the data shifted by the time interval of Ti. 

5. The interpolation circuit according to claim 4, 
characterized in that at least one of said fifth and sixth 
convolution operation unit comprises data holding unit for 
holding n pieces of data output from said oversampling 
operation unit while shifting, and addition unit for adding 
n pieces of data held in said data holding unit. 

6 . The interpolation circuit according to any one of claims 
1 through 5, characterized by comprising data appending unit 



for appending data having the symmetrical values proportional 
to the input data before and after said input data in the former 
stage of said oversampling operation unit. 
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ABSTRACT 



An interpolation circuit capable of performing 
interpolation operation with a simple constitution. A 
16 -times oversampling from discrete data is performed by D 
frip-frops 4,5. A first convolution operation is performed 
by D frip-frops 4 through 11 and an adder 12, on the result 
of which a second convolution operation is performed by D 
frip-frops 13 through 20 and an adder 21. Data interpolated 
along a quadratic function curve interpolating the discrete 
data is obtained from the adder 21. 
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This international preliminary examination report has been prepared by this International Preliminary Examining Authority 
and is transmitted to the applicant according to Article 36. 



This REPORT consists of a total of 



. sheets, including this cover sheet. 



| | This re P° rt is also accompanied by ANNEXES, i.e., sheets of the description, claims and/or drawings which have 
been amended and are the basis for this report and/or sheets containing rectifications made before this Authority (see 
Rule 70.16 and Section 607 of the Administrative Instructions under the PCT). 



These annexes consist of a total of 



sheets. 



This report contains indications relating to the following items: 
Basis of the report 
Priority 

Non-establishment of opinion with regard to novelty, inventive step and industrial applicability 
Lack of unity of invention 

Reasoned statement under Article 35(2) with regard to noveltv. inventive step or industrial applicability- 
citations and explanations supporting such statement " >f 

Certain documents cited 

Certain defects in the international application 
Certain observations on the international application 



I 


12<J 


II 


□ 


III 


□ 


IV 


□ 


V 




VI 


□ 


VII 


□ 


VIII 





Date of submission of the demand 
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Form PCT/IPEA/409 (cover sheet) (July 1998) 



INTERNATIONAL PRELIMINARY EXAMINATION REPORT 



International application No. 

PCT/JPOO/03040 



I. Basis of the report 



1. With regard to the elements of the international application:* 
the international application as originally filed 
the description: 

pages ^ 

pages 
pages 



□ 



. as originally filed 

. filed with the demand 



, filed with the letter of 



□ 



the claims: 
pages 

pages 

pages 
pages 



, as originally filed 

f as amended (together with any statement under Article 1 9 
, filed with the demand 



□ 



, filed with the letter of 



the drawings: 
pages 
pages 
pages 



, as originally filed 

, filed with the demand 



filed with the letter of 



the sequence listing part of the description: 
pages 
pages 
pages 



, as originally filed 

, filed with the demand 



, filed with the letter of 



2. With regard to the language, all the elements marked above were available or furnished to this Authority in the language in which 
the international application was filed, unless otherwise indicated under this item. ' 

These elements were available or furnished to this Authority in the following language which is- 

the language of a translation furnished for the purposes of international search (under Rule 23.1(b)). 
| 1 the language of publication of the international application (under Rule 48 3(b)) 

□ the language of the translation furnished for the purposes of international preliminary examination (under Rule 55 2 and/ 
or 55.3). 

With regard to any nucleotide and/or amino acid sequence disclosed in the international application, the international 
preliminary examination was carried out on the basis of the sequence listing- 

□ contained in the international application in written form. 

ED filed together with the international application in computer readable form 

□ furnished subsequently to this Authority in written form 

□ furnished subsequently to this Authority in computer readable form. 

ED The statement that the subsequently furnished written sequence listing does not go beyond the disclosure in the 
international application as filed has been furnished. 

EH The statement that the information recorded in computer readable form is identical to the written sequence listing has 
been furnished. 

The amendments have resulted in the cancellation of 

□ 

the description, pages 

□ the claims, Nos. 

□ 

the drawings, sheets/fig 



| | Thls re P° n has been established as if (some of) the amendments had not been made, since they have been considered to go 
beyond the disclosure as filed, as indicated in the Supplemental Box (Rule 70.2(c)).** 

* Replacement sheets which have been furnished to the receiving Office in response to an invitation under Article 1 4 are referred to 
l and70?7)° rt ^ ° nginally f iled " and are not annexed to this re P°rt since they do not contain amendments (Rule 70.16 

** Any replacement sheet containing such amendments mast be referred to under item 1 and annexed to this report 
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V. Reasoned statement under Article 35(2) with regard to novelty, inventive step or industrial applicability; 
citations and explanations supporting such statement 



Statement 
Novelty (N) 

Inventive step (IS) 
Industrial applicability (IA) 



Claims 
Claims 

Claims 
Claims 

Claims 
Claims 



1, 2, 4-6 



1, 2, 4-6 



1-6 



YES 
NO 
YES 
NO 

YES 
NO 



Citations and explanations 

Document 1: Qi Wang et al . , "Circuit Design of a D/A 
Converter Using Spline Functions", Signal 
Processing, March 1989, Vol. 16, No. 3, pp. 
279 to 288 



Document 2 



Document 3 



Akira Sakurai, "Supurain Kansu Nyumon", 
Tokyo, Tokyo Denki Daigaku Shuppan Kyoku, 
1986, page 7, line 21 to page 8, line 22; 
page 61, line 3 to page 103, line 1 

JP, 8-330957, A (Kenwood Corp.), December 13, 
1996 (13.12.96), paragraph [0012]; Fig. 2 to 
4 



Claims 1, 2, and 4 to 6 

Document 1 discloses a circuit for determining 
interpolation functions using calculations based on a 
number of finite observation values when individual values 
of interpolation functions which interpolate an infinite 
number of observation values are being determined. 

Moreover, Document 2 discloses the feature of 
determining spline functions by making the step functions 
indefinite integrals multiple times and of determining the 
spline function which passes through the data points by 
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using the aforementioned spline functions. 

However, an interpolation circuit which repeats a 
convolution operation more than twice with respect to 
zero-order-hold input data and generates interpolation 
data along a quadratic function curve that passes through 
the integral multiple of the input data value is neither 
disclosed nor suggested in any of the documents cited in 
the international search report and Document 3 newly cited 
in this written opinion. 

Claim 3 

Document 3 discloses a D/A conversion device provided 
with a sampling circuit (23) which resamples the digital 
data DATAp and outputs digital data DATAq in which the 
sampling period is an integral multiple, and an 
interpolation means which delays the data by latching in 
sequence the obtained DATAq to a plurality of latch 
circuits and forms linear interpolated data by adding this 
delayed data to the DATAq value. Since the invention set 
forth in Claim 3 forms a part of the D/A conversion device 
disclosed in the above-mentioned Document 3, Claim 3 lacks 
novelty . 
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VIIL Certain observations on the international application 

The following observations on the clarity of the claims, description, and drawings or on the question whether the claims are fully 
supported by the description, are made: 

The "convolution operation" disclosed in Claims 1 to 3 

is represented by the data generated as the result of the 

calculation with respect to the input data, but it is 
unclear what kind of operation this is. 

Claim 1 discloses the feature wherein the convolution 
operation is repeated twice or more and interpolation data 
is generated along a quadratic function curve. However, 
there are no disclosures in the description relating to 
the case when the convolution operation is repeated three 
times or more. Therefore, the description does not fully 
support the aforementioned feature set forth in Claim 1 
wherein the convolution operation is repeated three times 
or more. (Data is only generated along a quadratic curve 
when the convolution operation is repeated twice and when 
the convolution operation is repeated three times or more, 
the data should be generated along a dimension curve that 
is cubic or greater.) 

Claims 2 and 3 use the term "approximately" and, 
consequently, Claims 2 and 3 are unclear. 

Claim 2 discloses the feature wherein a second 
convolution operation is carried out with respect to zero- 
order-hold input data to determine a plurality of second 
data that forms a symmetrical trapezium. However, in the 
event that multiple sets of input data with different 
values are continuously inputted (see Fig. 11 (1)), the 
continuous input data is added to one another and 
therefore, it would not be possible to form a symmetrical 
trapezium (see Fig. 11 (5)). Consequently, it is not 
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VIIL Certain observations on the international application 

possible to equate the feature disclosed in Claim 2 
concerning the continuous input of multiple sets of input 
data with different values with the disclosures in the 
description and the figures. 

The same argument can be applied to the feature 
disclosed in Claim 3 wherein "the base of the isosceles 
triangle is approximately double the width of the 
aforementioned input data". 
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